Carbon Nanostructures as an Electromechanical Bicontinuum 
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A two-field model provides an unifying framework for elasticity, lattice dynamics and electrome- 
chanical coupling in graphene and carbon nanotubes, describes optical phonons, nontrivial acoustic 
branches, strain-induced gap opening, gap-induced phonon softening, doping-induced deformations, 
and even the hexagonal graphenic Brillouin zone, and thus explains and extends a previously dis- 
parate accumulation of analytical and computational results. 
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Vibrations in carbon nanostructures such as tubes, 
fullerenes, or graphene sheets [3, [3,0] have a ubiquitous 
influence on electronic, optical and thermal response: 
scattering from optical phonons limits charge transport 
in otherwise ballistic nanotube conductors 0, [B| ; twist 
deformations gap metallic tubes d, 0] ; ballistic phonons 
transport heat in nanotubes with great efficiency [H, 0, 
Hoj : resonant Raman spectroscopy can unambiguo usly 
identify a tube's wrapping indices (n,m) [IJ, 3, 
electron-phonon interactions may ultimately limit the 
electrical performance of graphene fltj| . Computa- 
tionally intensive atomistic models of lattice dynamics 
often lack simplified model descriptions that can facili- 
tate insight, yet traditional analytical continuum mod- 
els US Ell El, wnne ver y useful an d important, cannot 
describe atomistic phenomena without phenomenologi- 
cal extensions [IH, 

El ED]. Altho ugh continuum models 
are restricted to long-wavelength physics, they have been 
used to describe atomic-scale phenomena in bulk binary 
compounds by incorporating a separate continuum field 
for each sublattice [23j: in graphene, two fields are neces- 
sary. Here we present an analytical "bicontinuum" model 
that represents the full atomistic detail of the graphenic 
lattice, including optical modes, nonlinear dispersion of 
in-plane phonons, electromechanical effects and even the 
hexagonal graphenic Brillouin zone, a construct generally 
held to be exclusively atomistic. 

Graphene decomposes into the two triangular sublat- 
tices of Fig. [1] We describe in-plane deformations of the 
sublattices via two fields, ^(x), v 1 (x), i = 1, 2, and their 
strain tensors u 1 ^ = d^u^ and w y = d^v 1 ^ . The density 
of clastic energy contains direct and cross terms: 



V[u, v] = d[u] + d[v] + c[u, v] 



(1) 



Six-fold symmetry of the sublattices implies isotropy of 
the direct terms 12411: 
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d[u] — p! u^Uij + — u'„t( 
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FIG. 1: The two sublattices (circles and squares) of graphene 
and the three unit vectors e' 1 ' used in the text. <f>, z are 
cylindrical coordinates of a tube, while ^ = n/6 — 8 C with 6 C 
the chiral angle. Also, anisotropic (u xx = u xy = 0, u yy = 2^, 



q x — £j), shear (i 



= 0, u xy = rj, q y — —lr\) strains. 



c[u,w] = 2 Vij + Xu]Vj 
+ a(u — v) 2 

- (3e ijk {u^ +««) (u k -v k ) 



(3) 



The tensor e^fe, which is invariant under C$ v , can be 
represented by the three unit vectors {e™} of Fig. Q] 
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Symmetry dictates the form of the cross term 



Only the last term in Eq. [5] is not invariant under gen- 
eral rotation. (In nanotubes, it depends on the helical 
angle 9 C : = -e^ zz = -sin(36l c ), e zzz = -e^ z = 

— cos(3# c ), where 4>, z are defined in Fig. Q]). This elastic 
energy density, the lowest-order approximation in both 
derivatives and fields, contains six parameters: \j! and A', 
being confined to one sublattice, describe next-neighbor 
interactions; the cross terms \x and A describe nearest- 
neighbor interaction; a describes the stiffness against rel- 
ative shifts of the sublattices; f3 determines the strength 
of rotational symmetry breaking and so carries the point 
group symmetry of graphene. These parameters are nor- 
malized to the sublattice surface density cr s , so that the 
elastic energy is W — J <r s V d 2 x. 

Taking \a s (u 2 + i> 2 ) as the surface density of kinetic 
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energy, the equations of motion read 



d i u% (v) ' 
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with the sublattice 2-D stress tensors 

= 2fi' u ij +X'S ij ul + 2fi v lj + A S lj 



'(») 



= 2// u y - 
-/3e« (t 



A' <5^£ 



2/z 



(5) 



(6) 



As expected, a determines the frequency of two degener- 
ate k = optical modes: wr 2 = 4 a. 

First, we briefly show that the usual macroscopic elas- 
tic energy of graphene and its Lame coefficients can be 
obtained from V by considering a static, uniform solu- 
tion of Eqs.[5]with identical deformations on both lattices 
with an internal displacement 2q % = u l — v l : 



2f = le, 
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i-.lm 
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where I = (3 /a is a characteristic length. Anisotropic 



(27 



') and shear (77 = u* 3 ') strains produce in- 



ternal displacements q x — £j and q y = —in (Fig.Q}. The 
elastic energy for uniform deformations W u = J V u a g d 2 x 
then simplifies to 



V u [u,q] 



where a a 



m h — 

a 



u,, + - ( Xr- ^- ) ' 



2a<z 2 -2/3ey fe uV, 
2cr s = 2.26 g cm" 



(8) 



is the surface density 

of graphene, fi R = fx + fj,' - -fp> = A + A' + 
the measurable Lame coefficients |24[ . Macroscopic prob- 
lems do not distinguish between the two sublatices; elim- 
inating q % in Eq. [8] through Eqs. H] and [7] we obtain 
the familiar, isotropic, macroscopic energy for graphene, 
V u = /'/,•'/'■''/,, + Xru\u 3 J2. In the long wavelength 
limit Eqs. [5] returns the familiar longitudinal and trans- 
verse speeds of sound in terms of the Lame coefficients: 



Xr, v\ 



The out-of-plane displacements u±(x) and vj_(x) do 
not couple with the in-plane u l , v l in the harmonic limit: 
invariance under simultaneous sign change of u± and v± 
prevents it, for flat sheets. Introducing 2p±(x) = u±(x) + 
vj_(x) and 2q±(x) = ux{x) —vx(x), Vx must be invariant 
under px — > px + L(x), L(x) a linear function in the 
plane, and thus, can contain only second (and higher) 
derivatives in p± . Symmetry dictates (cf . Appendix) 

V x = laxql- Aa' ± d t qxd l qx+ Wxe ijk d k q ± d ij px 
+ 2(i+d ijP xd»px + X+dtpxdlpx 
- 2[j,±d ij qxd ij qx- X^diqxdjqx- (9) 

The frequency of the k = out-of-plane optical mode is 
2^/ax, and the out-of-plane acoustic branch is quadratic 
at small wave- vector, as expected. 
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FIG. 2: Bicontinuum phonons compared to EELS data (dia- 
monds [111 and squares [13]), fitting either to the entire Bril- 
louin zone (top) or just around F along F — » M. 



The bicontinuum phonons are much more richly struc- 
tured than in a traditional continuum model: they in- 
clude all the optical branches, show nonlinear dispersion 
at large wavevector, and even display the main features 
of the Brillouin zone, all without sacrificing the advan- 
tages of a continuum framework. Plane- wave solutions of 
Eqs. [5] returns an analytically solvable fourth-order secu- 
lar equation in w(k), yielding two acoustic and two opti- 
cal branches. The longitudinal branches cross at the ver- 
tices of a hexagon. Since the two-field elastic energy den- 
sity respects the point group symmetry of the graphene 
lattice, this hexagon is oriented just as the graphene Bril- 
louin zone; althou gh t he model, unlike in the envelope 
function approach [25|, has no built-in length scale, the 
elastic parameters can be constrained so that the crossing 
point coincides with the K point of graphene. A simi- 
lar argument holds for the out-of-plane modes: strikingly 
one can construct the correct Brillouin zone within a con- 
tinuum model. Fig. [5] shows the bicontinuum phonons fit 
to electron-energy- loss spectroscopy (EELS) data 26|, 27 1 
for parameters fitted either to the full Brillouin zone or 
just around T (28l |. 

The bicontinuum provides a unified framework for nan- 
otube mechanics which can describe all current compu- 
tational results on the coupling of nanotubc phonons 
to static structural distortions, to each other (e.g. 
breathing-to- Raman or longitudinal-to-transverse modes 
in helical tubes) and to the tube electronic structure. 
In a cylindrical geometry with coordinates {r, 0, z}, a 
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coupling between the tangential displacements u l , and 
the radial u r — u± appears in V of Eq. [1] via — 
) jr (and similarly for v); this accounts for 
the emergence of the Radial Breathing Mode (RBM) [2i| . 
We consider uniform solutions: u = u a e~ luJt , v — v e~ lu}t . 
The tube's helicity can be subsumed into new axes {£, £} 
(£ = 0cos30 c + zsin30 c ( = — 0sin30 c + zcos30 c ) ro- 
tated by an angle 30 c with respect to the base of the 
tube. In terms of p, q we obtain p^,p^ = and 



q< (u 2 



4a) + 2§p r = 

2 vl+P 2 /c*\ ^ 2 £<^ 



4a) = 
- 4a j_ 





2u-2u'+A-A' 




(10) 



Unlike standard elasticity [17j , which cannot describe op- 
tical modes, or standard atomistic descriptions, which 
cannot be solved analytically, the two-field continuum 
model enables an exact analytical solution for the cou- 
pling between the RBM and the graphite-like optical 
mode through the first two of Eqs. in (fTU)) : the RBM 
induces a shear in the sublattices, = tp"? = u r /r, 
which couples with the internal displacement through 
/3, and vice versa. Thus, the RBM is not purely ra- 
dial, but has a longitudinal component q z B ~ rt cos30 r: , 
as previously seen in a numerical calculation[30]. Ex- 
pansion of the RBM frequency in powers of l/r re- 
veals a correction to the the standard continuum re- 



sult vl/t [lTf: u>b 
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The 

/4a, 



graphite- like optical modes of chiral tubes are tu^ 

u>(/uj£ = 1 + | (£) + O , also of mixed longitudi- 
nal/transverse character except for armchair and zig-zag 
nanotubes, while the out-of-plane optical mode u± = 



4a i - 



2/J-2//+A-A' 



1/2 



is purely radial. A density 
functional theory calculation of the breathing mode [3l[ 
reports different frequencies with (ujb) and without (ujb) 
coupling to optical modes. We predict r 2 (ui B — w|) — > 
(3 2 /a as r — > oo: using ref [3l[ data for Qb, <^b we ob- 
tain I = PI a = 0.25 A (0.27 A) for non metallic zig-zag 
(armchair) tubes, in good agreement with the parameters 
from our fit to the graphene phonons [28| . 

The bicontinuum can also describe electron-lattice cou- 
pling to both acoustic and optical modes, by incorporat- 
ing a tight-binding model whose nearest neighbor hop- 
ping integrals t^ x \ t^ 2 \ are modulated by the in-plane 
elastic deformations: 



(11) 



where e is the inter-atomic distance and r a parameter 
to be determined [32j. For example, lattice deformations 
open gaps in metallic tubes, and these gaps in turn affect 
vibrational frequencies. If e c , e v are the conduction and 



valence bands, we have to nearest neighbors 



e c {kf-e v {kf 
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tW+2^t('¥ m )cos(fc-a (n) ), (12) 



where = — e^ m \ n(l,m) is cyclic in {1,2,3} (e.g. 

;( 2 )) and {eW} connects nearest neighbors. 



,0) 



= (1) 



From Eqs. 1111121 we find the band gap opened by strain 
in a metallic nanotube to be 



A 



(3t)' 



+ \ (ztq 1 ) 2 + - e ijk ^ h u ij 4> h q h - - (',.,/.//'-'o' v j . 

(13) 

In the second line of equation (|13[) the symmetry of the 
honeycomb lattice is broken by the unit vectors <p l , z l of 
the cylindrical coordinates. In terms of 27' = — u zz , 
if = u^ z , q z , equation (fT3"|) reads 

A = 3r|,f /e + 7'cos(36> c ) + r/sin(30 c )| , (14) 

which corrects and extends a well known previous result 
within a one-field continuum model p} that neglected the 
inner displacement (i.e. q l = 0). 

Opening bandgaps in metallic nanotubes causes several 
shifts in observed quantities. The term proportional to 
q 2 in Eq. (fT3|) show that longitudinal optical modes open 
a bandgap in metallic tubes of any helicity; the elastic en- 
ergy lowers by a term proportional to the square of the 
bandgap, leading to a the softening of longitudinal optical 
frequency in metallic nanotubes, as revealed by a recent 
DFT study [33j . Eq. (fl"3"| predicts also a softening of the 
RBM in metallic nanotubes = -A cos 2 (30 c ), high- 
est for zig-zag tubes as seen in DFT [3l| , and relates it to 
the optical softening, with A = (1 — 1/ e)ui opt Suj pte 2 / iv^, 
tOopt the graphite-like optical mode, and 8<jj op t its soften- 
ing in metallic tubes (A ~ 2%). Other shifts can be 
predicted: the speed of sound for the twist mode softens 
by = -^-Asm 2 (30 c ), or ~ 2.2% in armchair tubes. 

Doping-induced structural deformations can also be 
studied by minimizing the total energy (elastic plus 
doped electrons). Subtle phenomena absent in other 
models [22| can be accessed within the bicontinuum 
framework. Going to next-nearest-neighbor in the hop- 
ping integrals (dt^ = — n a\ afpv %: > 32]), we find that 
at first order in both a/r and the number of dopant elec- 
trons per atom p e , semiconducting (n, 0) nanotubes show 
doping-induced changes in tube length (dL/L = u zz ) and 
axial bond-length (db ax — eit 22 — q z ): 




_PcT 



± 



i 1 



?\ r 3ri 2/j.n+Xn 
') f 2t hr+Xb 



(15) 



where mc is the mass of the carbon atom. The sign is 
positive (negative) for r = n mod 3 = 2 (n mod 3 = 1). 
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Recent DFT results (34( indeed show shrinking or stretch- 
ing of b ax for n = 16,13 or n = 14,11 tubes respec- 
tively, as predicted by Eq. [T5j In DFT, the overall tube 
lengthens in the second case (n — 14, 11), again in ac- 
cord with the bicontinuum; the lengthening found for 
r = 2, is less than for r = 1, perhaps a consequence of 
the change in sign in Eqs. 1151 Finally the shrinking of 
the axial bond determines an up-shift in the longitudinal 
graphite-like optical mode and might explain recent Ra- 
man results that point toward anomalous bond contrac- 



tion under doping in semiconducting nanotubes [351 . 

In summary, a symmetrized two-field continuum model 
of graphene and carbon nanotubes provides the first uni- 
fied analytical treatment for a wide range of vibrational 
and electromechanical phenomena including nonlinear 
dispersion of in-plane phonons, zone-edge degeneracies 
and optical modes. A full range of vibrational-electronic- 
mechanical couplings, which were absent from previous 
continuum models or happened upon in an ad hoc fashion 
in computational work, can now be understood within a 
single unified analytical framework. Extending the for- 
malism to include higher-order effects arising from curva- 
ture or metallic character (i.e. symmetry breaking terms 
containing (jf , z l , as in Eq. I13|) . anharmonicity (terms 
higher order in u lJ , v lJ ), or long-distance interactions 
(higher partial derivatives) is straightforward. An ex- 
tension to boron nitride nanotubes, with different coef- 
ficients for each sublattice in the direct terms of Eq. [2l 
might prove useful to study their piezoelectricity. 

Appendix: Derivation of Eq. \3\ 

The term c[u, v] must be invariant under the combina- 
tion of 2n/6 rotations and the exchange of fields u <-> v. 
Adding reflection through the x axis (Fig. [J) then im- 
plies Cs v invariance. There is also a field translation 
invariance: u(x) —* u{x) + p, v(x) — > v(x) + p. The 
objects u\ v 3 , u lJ , and v lJ can be combined pairwise 
only into tensors of rank two, three and four; thus c[u, v] 
decomposes into three parts. The first part has terms 



u v- 



symmetry then implies the form a(u 



like 

with a > to ensure an energy minimum. The sec- 
ond part has terms like u % ^v kl ; the only admissible form 
is 2/itt lJ Uy- + \u\Vj. The third part contains only rank 
three terms such as u^v k contracted with a C^ v invari- 
ant tensor e^fc, giving eijkU % ^v k . By requiring invariance 
under 27r/6 rotations conjugated with sublattice switch- 
ing, and also the field translation invariance, we obtain 
the form e i3 kU lj (u k — v k ) + e* jk v lj (v k — u k ), where the 
star means a 2-7r/6 rotation. Since C^ v invariance implies 



-ijk 

[1] 
[2] 



e-ijk we finally obtain the third row of Eq. [3l 
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